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ABSTRACT 

Virtually all detailed gas-phase models of the chemistry of dense interstel- 
lar clouds exclude polycyclic aromatic hydrocarbons (PAH's). This omission 
is unfortunate because from the few studies that have been done on the sub- 
ject, it is known that the inclusion of PAH's can affect the gas-phase chemistry 
strongly. We have added PAH's to our network to determine the role they play 
in the chemistry of cold dense cores. Initially, only the chemistry of neutral 
and negatively-charged PAH species was considered, since it was assumed that 
positively-charged PAH's are of little importance. Subsequently, this assumption 
was checked and confirmed. In the models presented here, we include radiative 
attachment to form PAH~, mutual neutralization between PAH anions and small 
positively-charged ions, and photodetachment. We also test the sensitivity of our 
results to changes in the size and abundance of the PAH's. Our results confirm 
that the inclusion of PAH's changes many of the calculated abundances of smaller 
species considerably. In TMC-1, the general agreement with observations is sig- 
nificantly improved contrary to L134N. This may indicate a difference in PAH 
properties between the two regions. With the inclusion of PAH's in dense cloud 
chemistry, high-metal elemental abundances give a satisfactory agreement with 
observations. As a result, we do not need to decrease the observed elemental 
abundances of all metals and we do not need to vary the elemental C/0 ratio in 
order to produce large abundances of carbon species in TMC-1 (CP). 



Subject headings: Astrochemistry - Molecular processes - ISM: abundances - 
ISM: molecules - ISM: individual objects: L134N,TMC-1 (CP) 
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Introduction 



Polycyclic Aromatic Hydrocarbons (PAH's) are liydrocarbon compounds witli multi- 
ple six-membered aromatic rings. Such species are thought to be a c ommon comp onent of 
the interstellar medium in the Milky Way and in external galaxies ( Drainell2003l ). PAH's 
are typically observed in infra-red emission bands following excitation by either ultra-violet 
or visible radiation in assorted regions of the interstellar medium exposed to radiatio n 
( Leger fc D'hendecourt 1985 : Leger et "al1ll989 : AUamandola et al. 1985 : Salama et al. 1999). 
Despit e all attempts, no specific types of PAH's responsible for these features have been iden- 



tified ( Mulas et al. 



20061). and the fo rm and intensity of the spectra vary strongly from one 



region to another (jPraine fc Lill2007l ). probably reflecting a v ariety of sizes a nd abundances. 
The origin of PAH's also remains an unanswered question (IWilliamsl l2003[). Most beli eve 
that PAH's are formed in the atm ospheres of carbon stars (lAUain et al.l 119971 : ICaijl2002l ) or 
in super novae (ITappe et al.ll2006l ). 



Various observations show that the abundance of small free-fiying P AH's have a ten- 



dency to decrease in dense clouds where UV p hotons do not penetrate (jBoulanger et al. 
199ol : [Bernard et al.lEggsl : iRapacioh et al.ll2005h . iRapacioh et~al1 J2005I . booeh proposed that 



PAH's form larger PAH's or PAH-clusters in such dense regions. As a consequence, free-fiying 
PAH's observed at the border of illuminated clouds may be the result of photodissociation 
of these larger species. Although they determined a minimal size of 400 carbons for these 
species using mid- infrared ISOCAM maps of photodissociation regions (PDR's), this number 
is very uncertain. If free-flying PAH's are not present in the dense interstellar medium, an 
efficient process is needed to aggregate them. Aggregation is poorly understood, and the 
time in which this process occurs and the fraction of small PAH's remaining in dense clouds 
unexposed to intense radiation are undetermined. Currently, constraints on the PAH size 
and abundance as a function of the physical conditions in the dense interstellar medium 
remain weak. In this article, we focus not on the size of PAH's but on the chemical impact 
of their presence in the dense cold interstellar cores. Nevertheless, the sensitivity of our 
results to size will be examined. We follow past authors by lumping all PAH's into three 
groups of species: neutral, positively charged, and negatively charged PAH's, designated 
PAH, PAH"*", and PAH~ respectively. Of these three groupings, PAH+ species were found 
to be unimportant and were excluded from the detailed calculations. 

The chemical effect of PAH's at steady-state has been studied previously. In the diffuse 
interstellar medium, the inclusion of PAH's, among other effects, enhances the ne utralization 



of at o mic ions such as carb on, silicon and sulfur, reducing their abund ances (ILepp et al. 



f|l988h 



19881: Bakes fc Tielens 1998). In their study of dense cloud chemistry, Lepp fc Dalgarno 



concluded that for PAH's consisting of molecules with 30-50 carbon atoms, with 
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a fractional abundances greater than 10~®, PAH~ is the main source of negative charge, 
which considerably changes the abundances of some neu tral species, enhancing, for example , 
the abundances of carbon-bearing molecules. Finally, iFlower fc Pineau des Foretd (120031 ) 
studied the importance of PAH's for C-shock models. In this article, we report a time- 
dependent study of the role of PAH's in the gas-phase chemistry of cold dense clouds and 
compare our results with observations. One focus of our study will be how the presence of 
PAH's affects the choice of elemental abundances. 

The elemental abundances typically used for chemical mo dels of dense clouds are re- 
ferred to by the appellation "low-metal". Originally defined by iGraedel et al.l (Il982l ). these 
abun dances are ba sed on the Copernicus satellite observations of the ( Ophiuchi diffuse re- 
gion (jMortonlll974l ). But, in order to reproduce the observation of molecules in dense clouds 
better, the ( Oph abundances of all the elements but. He, C, N and O are decreased by 
two orders of magnitude. This assumption implies that a strong depletion of the heavier 
elements occurs between the diffuse and dense cloud phases (IGraedel et al.lll982l ). Although 
the details of this depletion are not understood, the depleted abundances are necessary to 
reproduce the abundances of molecular ions and, to a lesser extent, carbon-bearing species. 
Since the presence of PAH's in dense clouds should modify the ionization fraction, it may 
help in solving the problem of the choice of elemental abundances. 

This article is organized as follows. In Sect. [2], we describe the parameters of the PAH's 
used for the analysis and the chemical processes involving them. In Sect. [31 other details of 
the chemical model are given. The influence of PAH's on the abundances of smaller species 
is reported in the next section, while the sensitivity of the results to the PAH parameters is 
discussed in Sect. [51 In Sect. 6, we compare observations in molecular clouds with grids of 
models in which we vary several parameters, including PAH size and abundance. The paper 
ends with our conclusions. 



2. PAH parameters and chemistry 



Draine fc Lil (120071 ) reproduced the extinction curves of various diffuse regions of the 



Milky Way from Spitzer observations and constrained the size distributio n of carbonaceous 
grains inc luding PAH's down to 3 A in size (see Fig. 6 of lDraine fc Lill2007l ). Using the model 
Jm = 1 of lDraine fc Lil (120071 ) . which has the largest abundance of small PAH's, we integrated 
over a radius range of 3.8 to 6.0J^. The integration yields a mean radius of 4 A (Nc = 30) 



^The dependence of the number of ca rbon atoms N g on radius a for PAH's is not well-known s ince i t 
depends on the structure of the molecule (Omont 19861 ). We used the relation given by Draine fc Li ( 2007 ): 
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and an integrated fractional abundance with respect to rin of 3.07 x 10~^. Multiplying by 
30 leads to a fractional abundance of C atoms in PAH's of 9.1 x 10^^, or roughly 3% of the 
elemental carbon abundance. 

Because PAH's may not be free in dense clouds, we dec i ded t o also consider PAH-clusters 
of radius 10 and 100 A. With the formula of IPraine fc Lil ( 20071 ). these limits correspond to 
PAH's with 468 and 4.68 x 10^ carbon atoms, leading to fractional abundances of 2 x 10~^ 
and 2 x 10~^^, respectively, if one assumes (a) the same atomic carbon abundance in PAH's 
as for the smaller species, and (b) that weakly-bound clusters can also be approximated by 
the same relationship between number of carbon atoms and size. 



Bakes fc TielensI ( 1l998l ) and lDartois fc D'HendecourtI (119971 ) f ound that the abundance 



of PAH"*" decreases strongly in clouds with an A„ larger than 2 and iLepp fc Dalgarnd (119881 ) 
found that the abundance of PAH"*" is smaller than the abundance of the other forms of 
PAH's by about one order of magnitude in their calculations. For this reason, in our study of 
dense cloud chemistry, we did not originally include cations but only neutral and negatively- 
charged PAH's. Prodded by the referee, however, we looked more carefully at the role of 
PAH"*", but concluded again that their exclusion is justified. The reconsideration will be 
discussed in t he next section. We st art t he chemistry with neutra l PAH 's only, based on 
the results of iBakes fc TielensI (Il998[ ) and lOartois fc D'HendecourtI Jl997h . Neutral PAH's 
become negatively charged very rapidly and our results after a very small time interval do 
not depend on whether the chemistry commences with PAH or PAH~. 



2.1. Chemical processes involving PAH's 



The main problem of including PAH's in chemical models is that few experimental mea- 
surements of reaction rates have been performed on th ese molecule s , and these few tend to 
involve rather small cationic PAH's only. For example, ISnow et al.l (119981 ) studied the rates 
of the reactions between a variety of positively-charg ed PAH's wi t h six, ten, and sixteen car- 
bon atoms and the neutrals H, O, N, and H2, while iBetts et al.l (120061 ) studied reactions of 
coronene cation (C24H^2) with the same reactants. Although association reactions appear to 
be dominant, most the experiments were done at sufficiently high densities that the mech- 
anism, three-body association, is not what occurs in the low-density interstellar medium. 
Moreover, cationic PAH's are not important in our study. 



Nc — 460(a/10A)'^. The radius range here then corresponds to PAH's with between 25 and 100 carbon 
atoms. 
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The processes in our model include electron attachment, recombination of negative PAH 
ions with smaller positive ions, and photodetachment of negative PAH ions by UV radiation. 
We consider these in turn. 



Negative PAH ions are formed by radiative electron attachment : 



PAH 



PAH~ + hu. 



(1) 



The r ate c oefficient for this reaction is thought to depend on the size of the PAH. lOmont 
(119861 ) and lAUamandola et al.l (Il989l ) proposed two different relations: ke = 10^'^ (Nc)^/'^ 
and fee = 1.2 X 10~^(Ai'c)^/^ cm^ s~^. These two formulae lead to differences of about a fac- 



—1 

cm s 



tor of 2 for PAH's with Nc < 100. We arbitrarily decided to follow lie Page et all (120011 1 
and use the first formula. The rate coefficient is 1.3 x 10~^ cm^ s~^ for the attachment of 
electrons to neutral PAH's with 30 carbons. The results depend on this value only if ke 
is decreased by a factor of 100. In other words, an increase of this rate coefficient does 
not chan ge the results. Calc u lation s by Herbst (in preparation) based on the statistical ap- 
proach of iTerzieva fc HerbstI (120001 ) indicate that with typical electron affinities in the range 
of 0.5 - 1.0 eV, the large rate coefficients deduced by previous authors are valid for neutral 
PAH's with 30 carbon atoms, but not for t hose significantly smaller, in agreement with the 
experimental results of iMoustefaoui et al.l (119981 ). 



Negatively- charged PAH's can recombine with positive ions, designated X+. These 
recombination reaction s differ from dissociati ve recombination of positive ions with electrons 
in two important ways (iBates &: Herbstlll988l ). First, they are less exothermic because of the 
electron affinity of the PAH. Secondly, some of the reaction exothermicity can be deposited 
in internal energy of the PAH so that less is available to dissociate the positive ion upon 
neutralization. We assume that two resulting channels are possible for these recombinations, 
simple charge exchange: 

PAH" + X+ — > PAH + X, (2) 
and, if the positive ion contains hydrogen: 



PAH" + X^ 



PAH + X' + H, 



(3) 



where X' designates a neutral with one fewer H atom than the ion. The terms "neutraliza- 
tion" and "dissociative neutralization" are used for these two processes to distinguish them 
from recombination with electrons. We assume equal branching fractions for the two sets of 
products, if they exist. If the neutralization leads to a neutral species that does not exist 
or is not included in the model, only the second channel is considered. The products of the 
dissociation are very uncertain, and laboratory experiments are required. 
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The rate coefficients for these reactions have been computed by taking into account the 

Coul omb attraction between the charged species via the expression flSpitzedll94ll : IWeisheit fc UphamI 



1978|): 



k = HG? (Vq) ( 1 + 



where e is the electronic charge, is the Boltzmann constant, a is the radius of the PAH 
(or PAH cluster), and (vo) is the mean relative velocity between PAH~ and ions, which is 
given by the expression 



/8kBT 

TT/i ' 



(5) 



with T the temperature and /i the reduced mass, which is mainly determined by the mass 
of X"*". In expression (jl]), we have assumed that the radius of the PAH is the critical radius 
for which reaction can occur. The radius at which the second term equals the ffist term 
on the right-hand-side is approximately 1 /im at 10 K. This radius is much larger than the 
PAH radii considered here, so that 



2 X 10 ^ cm^ s^ 



aksf >> 1- The rate coefficient at 10 K is approximately 
for an ion of 30 amu and a PAH of 4 A. 



Relation (jl]) is based on the fact that, due to the delocalization of the negative charge on 
the PAH, the trajectory of the positive ion as it approaches the PAH will always be affected 
by this charge. Taking into account the complex shape and large number of possible PAH 
clusters in dense clouds, several negative charges may stick, and the attraction between such 
a system and a positive ion may depend on the localization of the electrons in the cluster. 
As a consequence, expression (jlj) may be more likely to be correct for small PAH's of 4A 
radius than for larger ones. For this reason, we consider the alternative possibility that rate 
coefficients for the reactions between PAH~ and X+ are given by the simple ballistic relation 



k = Tia^ (vq) . 



(6) 



For small PAH's, the rate coefficient varies as a whereas for large PAH's, if the process is 
ballistic, it varies as a^. For a 100 A PAH and an ion of 30 amu, the rate coefficient at 10 K 
is then approximately 2 x 10"^ cm'^ s~^, or 0.1 of the value for the smaller PAH's. 

In addition to these reactions, PAH~ is neutralized by photodetachment of electrons 
due to ultraviolet, visible, and even near infra-red photons: 



PAH" + hz/ 



PAH + e" 



(7) 



Following iLepp et al.l (Il988l ). we assume a cross section of 10 cm significantly above 



threshold, which is the electron affinity (see also iMalloci et al.l 120071 ). For a standard in- 
terstellar radiation field, this leads to a conservative estimate for the rate coefficient in the 
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absence of shielding of approximately 10~^ s~^, neglecting differences in electron affinities 
among PAH's. A crude estimate of the dependence on is given by the simple exponent 
exp{—A^). Even with an of 10, the rate is not negligible in dense clouds. Indeed, under 
most conditions, photodetachment is more rapid than neutralization, but it must be remem- 
bered also that the electrons produced in photodetachment then attach to neutral PAH. We 
do not consider photodetachment by cosmic ray-induced photons. If w e assume a rate coeffi- 



cient of 2 X 10"^ X ( for this process (IFlower fc Pineau des Foretsll2003l ). the photodetachment 



by direct UV photons is then more efficient by one order of magnitude even at an A^ of 10. 

In our original calculations, we did not include PAH+. In dense clouds, positively 
charged PAH's can be produced by secondary UV photons induced by cosmic-rays and 
by charge transfer with small positive ions X"*". The latter process can neutralize the X"*" 
abundance and so affect the small-molecule chemistry. Then PAH"*" can be destroyed by 
recombination with free electrons. As suggested by the refere e, we tested the importance 



of pos itively-charged PAH species using rate coefficients from iFlower fc Pineau des Forets 



( I2OO3I ) for PAH"*" production by photoionization and for recombination with electrons. For 



the charge-transfer reactio ns, we used ra te coefficients computed with the Langevin theory 



and the polarizability from lOmontI (119861 ). As an example, the Couloub-based rate coefficient 
for the recombination reaction between PAH~ and X+ is 7.5 x 10~^ cm^ s~^ for PAH's of 4 
A and an ion of 30 amu whereas the rate coefficient for the charge-transfer reaction between 
PAH and X^ for the same ion is 3.2 x 10^^ cm^ s~^. In the various models that we ran, we 
found that the abundance of PAH+ is at least two orders of magnitude smaller than the PAH 
and PAH~ abundances. Moreover, its inclusion does not affect any of the results discussed 
below. This conclusion holds even when the recombination rate coefficient is the ballistic 
one. 



3. Model parameters 

We utilize a pseudo-time-dependent gas-phase model, which computes the chemical 
evolution at a fixed temperature and density. A version of this model (NAHOON) is now 
available online at 

http : / /www . obs . u-bordeauxl . f r/radio/VWakelain/Valentine\y„20Wakelain/Downloads . html. 
The chemical database is a recent version of the OSU network: osu_01_2007 , which can be 
found online at http://www.physics.ohio-state.edu/~eric/research.html. The reac- 
tions involving PAH's will also be available at the same address. 

We use three different sets of elemental abundances (EAl, EA2 and EA3), which are 
listed in Table 1 with respect to the total nuclear hydrogen abundance along with the initial 
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choice of ionization. EAl represents the low-metal abundances discussed earlier. EA2 is 
based on the high-metal abundances observed in the ( Ophiuchi diffuse cloud, but is modified 
based on recent observations. The abundance of helium has not been observed in cold cores 
and we use a value of 0.09 for EA2 with respect to hydrog en, based on observations in the 
Orion Nebula (IBaldwin et al.l Il99ll : lOsterbrock et al.l Il992l ) . To account for the oxygen and 
carbon that are dep leted on grain mantles in the form of H2O and CO, we removed 2.4% 
of C and 13% of O (IShalabiea fc Greenberglll995l ). It is possible that the metals are more 
depleted in dense clouds than in diffuse regions. Thus, for EA3, we assumed that the silicon 
and magnesium are tota lly depleted on grains and we assum ed an abundance of 1.5 x 10^* 
for iron as suggested by iFlower fc Pineau des Foretd (120031 ) . The cor nplete depletion of S i 
is one explanation for the non-detection of SiO in dense clouds (see IZiurys et al.l Il989bl ). 
We also considered an additional depletion of a fact or of 10 for chlorin e as suggested by the 
observation of HCl in cold material around OMC-1 (IBlake et al.lll985l ). Note that in all the 
elemental abundances that we considered, the elemental C/0 ratio is 0.4. 

The species are assumed to be initially in atomic form except for hydrogen, which is 
assumed to be entirely in its molecular form. We consider typical physical conditions for 
dense clouds: a temperature of 10 K, an H2 density of 10^ cm~^, a cosmic-ray ionization rate 
( of 1.3 X 10~^^ s~^ and a visual extinction of 10. 



4. General effects of PAH's 

In this section, we present results for the three sets of elemental abundances (EAl, EA2 
and EA3) with and without PAH's. In the models with PAH's, referred to as EA/PAH, we 
assume their fractional abundance to be 3.07 x 10"'', our standard value for small PAH's. 
Variations in this value will be considered in Sect. [5] and [61 



4.1. Ionization fraction 

Figure [T] shows the ionization fraction (sum of the negatively or positively charged 
species) as a function of time for the three sets of elemental abundances, each with and 
without PAH's. The initial ionization is typical of a diffuse cloud since we start with ionized 
carbon in all models. With low-metal elemental abundances (EAl) and no PAH's, the 
ionization fraction stays high for only 10'^ yr and then decreases sharply, reaching a steady- 
state value of 3 X 10~® by 10^ yr. With both EA2 and EA3, and no PAH's, the ionization 
fraction decreases slowly, remaining larger than 10~^ for about 10^ yr and eventually reaching 
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a steady-state value of a few 10~^. When PAH's are added, the ionization fraction does not 
depend significantly on the elemental abundances. The main effect is to sharply decrease the 
earlier ionization fraction of the gas: within the first 300 yr, the ionization fraction goes from 
~ 10~^ to ~ 5 X 10~*, its steady-state value, which is only slightly larger than that achieved 
with the EAl model. Thus, the inclusion of PAH's reduces the ionization fraction for EA2 
and EA3 to a value close to that achieved with low-metal abundances and no PAH's. The 
steady-state ionization fraction for model EAl/PAH is larger by a factor of 2 compared with 
model EAl without PAH's. The ionization fraction of model EA2 /PAH at steady-state can 
be compared with the very similar value of 3.8 x 10^^ reported by lLepp &: Dalgarno Jl988h 



for a ratio of total PAH to hydrogen of 10 and high-metal abundances. 

Figure [2] depicts the fractional abundances of neutral and negatively- charged PAH's 
as well as the electron abundance as functions of time for model EA2/PAH. The other 
PAH models are not shown because these abundances do not depend significantly on the 
chosen elemental abundances. By 300 yr, the fractional abundances of electrons and PAH~ 
drop to 10~^ and ~ 7 x 10~^ respectively, so that PAH~ becomes the dominant negatively- 
charged species. The conversion of neutral to negatively- charged PAH's occurs via electron 
attachment, which is very efficient at early times because of the high electron abundance. 
But, once formed, PAH~ is efficiently destroyed by reactions with atomic ions, which are 
also very abundant at the beginning of the chemical evolution. Note that the recombination 
of atomic ions is much more efficient with PAH's than with electrons since the latter must 
occur radiatively. For instance, the neutralization of by PAH^ (with a = AA) is four 
orders of magnitude larger than by electrons at 10 K. Because of this efficient neutralization, 
the ionization fraction of the gas decreases much more quickly when PAH's are present (see 
Fig. [1]). We also find that the abundance of neutral PAH's exceeds that of the anionic 
PAH's b y a factor of 5-10 after 3 00 yr, which is also in agreement with the steady-state 



results of iLepp fc Dalgarnd (119881 ). 



4.2. Results for selected species 



In the absence of PAH's, the higher-metal elemental abundances EA2 and EA3 lead to a 
very different chemistry from EAl, and tend to worsen the agre ement between the m odeling 
results and the observations for selected classes of molecules (iGraedel et al.lll982l ). With 
models EA2/PAH and EA3/PAH, however, the chemistry is much closer to that computed 
with the low metal elemental abundances in the absence of PAH's. The major change in the 
chemistry is caused by the high abundance of PAH~, which neutralizes atomic positive ions 
efficiently and, upon reaction with positive molecular ions, tends to neutralize them without 
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significant dissociation, contrary to dissociative recombination reactions with electrons. In 
addition, the decrease of the ionization fraction at initial stages of the chemistry increases the 
abundance of ionic species such as H3" and HsO"*" because they are depleted more slowly. At 
later times, some of the C-, S-, N-bearing species are enhanced in the models EA2/PAH and 
EA3/PAH because of i) larger elemental abundances and ii) a smaller electron abundance. 
We now present the results for the most important species and for selected molecules that 
show significant differences among models. 



4.2.1. C, and CO 



In Figure [3l we show our results for C and C"*" as functions of time for the three elemental 
abundances with and without PAH's. As expected, the fractional abundance of ionized 
carbon is decreased most rapidly when the PAH's are included. The agr eement at steady- 
state f or the model EA2/PAH is in reasonable agreement with the result of iLepp fc Dalgarno 
(119881 ). A consequence of the C"*" neutralization by PAH~ is that atomic carbon becomes 
the main reservoir of carbon earlier than with the models without PAH's. At the age of 
dense clouds, between 10^ and 10^ yr, atomic carbon is more abundant in EA2 and EA3 
compared with EAl because of a higher elemental abundance of carbon. Our steady-state 
fractional abundance of neutral carbon (~ 3 x 1 0~^°) in model EA2/PAH is three orders of 
magnitude below that of the analogous model in lLepp fc Dalgarnd (119881 ). After early time, 
CO becomes the main repository of carbon in all models; the CO abundance (not shown in 
the figure) is insensitive to the model used. 



4.2.2. H+, HCO+, BiO+, OH, H2O and O2 

The rapid decrease of the electron abundance in the models with PAH's leads to strong 
increases in the abundances of molecular ions such as H^, HaO^ and HCO"*" and neutral 
species formed from them at the initial stage of the chemistry, as can be seen in Fig. |H 
At longer times, the ionic abundances show much less time dependence, with two general 
but not universal patterns: (i) they are more abundant in the models with PAH's than in 
the analogous models without them, and (ii) they are most abundant in models with EAl 
abundances and least abundant in models with EA2 abundances, with and without PAH's. 
General reasons for the trends include the higher electron abundances in models EA2 and 
EA3, the slower depletion rates by negatively charged species (mainly PAH~) in the PAH- 
containing models, and differences in abundances of important neutral species involved in 
depletion processes among the PAH-containing models. 
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The ion H.^ is a partial exception to the trends in that its abundance is not enhanced 
by PAH's at later times, mainly because it is destroyed principally by heavy neutral species 
such as CO and H2O rather than by charged species. The abundance in EA2 and EA3 
is smaller than in EAl because of dissociative recombination by more abundant electrons 
whereas it is smaller in EA2/PAH and EA3/PAH than in EAl/PAH because of reaction 
with H2O, which is more abundant in the high- metal models. The enhancement in the ion 
HCO^ for models with PAH's in the initial stages of the chemistry is even more dramatic 
than that of Hjj" because of its formation from abundant carbon atoms, by reactions leading 
first to CHj^. At later times, HCO"*" is formed by the reaction between H;^ and CO in all 
models. Without PAH's, HCO"*" is destroyed by electrons whereas it is destroyed by H2O in 
models with PAH's. 

The abundance of HsO"*" is also dramatically enhanced in the initial stages of the chem- 
istry in the presence of PAH's. Without PAH's, HsO^ is destroyed by electrons whereas it 
is destroyed by PAH" in EAl/PAH and EA3/PAH. In EA2/PAH, due to a large abundance 
of elemental silicon, HaO^ reacts with gas phase SiO, which is as abundant as 1.5 x 10~^ at 
10^ yr. Such a large abundance of SiO is in disagreement with observations in dense clouds 



fIZiurvs et al.lll989bh . 



One consequence of the large HsO"*" abundance, at initial stages of the chemistry, is that 
the abundances of OH, H2O and O2, which stem from the recombination and neutralization 
of HaO^, are also enhanced in PAH-containing models . The situation is somewhat more 
complex at later times with only water showing a much higher abundance in PAH-containing 
models. The predicted enhancement of H2O in PAH-containing models is a problem because 
the calculated abunda nces far exceed th e observational upper limits of ~ 3 x 10~^ in L134N 



and 7x 10"*^ in TMC-1 (jSnell et al.ll2000l ). The result for molecular oxygen is, on the contrary, 
not affected by PAH's at times after 10^ yr. Our theoretical O2 abundance is about 10"^ 
at 10^ yr using EA2/PAH or EA3/PAH, which compares favo rably with the obse rvational 
upper limits of 1.7 x 10"^ in L134N and 7.7 x 10"^ in TMC-1 dPaeani et al1l2003h . 



4-2.3. Assorted Hydrocarbons and Oxygen-containing Organic Molecules 

The computed abundances of the hydrocarbons C-C3H2 (cyclopropenylidene), C3H4 
(methyl acetylene), C6H2 (triacetylene) and CH3C4H (methyl diacetylene) are shown in 
Fig. [5] as functions of time. Without PAH's, the abundances peak at so-called early-times 
(10^ yr) for low- metal abundances (EAl), a well-known effect. For EA2 and EA3, on the 
other hand, the early-time peaks are weak to non-existent and the peak abundances very 
low, an effect caused by the high electron abundances, which disrupt the growth of molec- 
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ular ions. When PAH's are added, the abundances of the four hydrocarbons depicted are 
enhanced for all three models. Although the enhancement is significantly greater for EA2 
and EA3, the peak EAl values are still the largest. The enhancements are caused mainly 
by the lesser destructive effect of PAH~ compared with electrons when it reacts with in- 
termediate ions in the chain of synthetic reactions, although in individual cases secondary 
effects can be isolated. The major effect, which comprises both lower rates and larger prod- 
uct neutral channels, is stronger than the enhancement in destruction rates of hydrocarbons 
due to often larger abundances of molecular ions. That differences remain among the PAH 
models despite the fact that PAH~ has the same abundance shows, once again, that other 
processes are also important. The comp uted abundances o f hyd rocarbons at steady-state 
can be very different from the results of iLepp &: Dalgarnol (119881 ) . For example, the C,sH 2 



abundance is ~ 10 in model EA2/PAH whereas it is ~ 10 ^ in lLepp fc Dalgarnol (119881 ). 



This difference may be due to discrepancies in the network for carbon chemistry. 

As can be seen in Fig. |6l the abundances of the oxygen- containing organic neutral 
molecules HCOOH, CH2CO, CH3OH and CH3CHO in models containing PAH's are en- 
hanced at almost all times by at least one order of magnitude over the abundances in models 
without PAH's. Once again, the reason is generally the same as for the hydrocarbons: 
neutralization reactions between precursor ions with PAH~ are not as destructive as disso- 
ciative recombination reactions with electrons. Even if a species, such as ketene, is formed 
by a neutral-neutral reaction, its precursor is likely formed at least partially through ion- 
molecule channels. An as example, consider the case of methanol. Methanol (CH3OH) 
is produced through the ion CH30H^, which itself is formed via the radiative association 
between CB.^ and water. In models without PAH's, the protonated ion can recombine dis- 
sociatively w ith electrons to form C H3OH + H, although this channel is now known to be 
a minor one (IGeppert fc et al.l 120061 ). With PAH~, on the other hand, we assume that this 
channel is the only one, since the neutral form CH3OH2 does not exist and we only allow one 
hydrogen to be ejected from the species (see section [^7T]) . Another reason for the enhanced 
abundance of methanol lies in the enhanced abundance of water. We note that with PAH's, 
the peak predicted fractional abundance of m ethanol is still 1-2 o rders of magnitude below 
the typical value in cold dense clouds of 10"^ ( IGarrod et al.l 120071 ). 



Unlike the situation with the hydrocarbons, the abundances in models EA2/PAH and 
EA3/PAH are not significantly less and sometimes even greater than those in model EAl/PAH. 
The small differences suggest that PAH~ plays an even more dominant role than in the case 
of the hydrocarbons. In the non-PAH models, although the EAl abundances tend to be 
greater than those in EA2 and EA3, the effect is smaller than for hydrocarbons except for 
acet aldehyde. 
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4- 2.4- Cyanopolyynes 



Cyanopolyynes are very abundant towards the famous cyanopolyyne peak (CP) of the 
cold core TMCl, yet their abundances are un der produced by gas-phas e models unless a high 
C/0 elemental abundance ratio is assumed (ITerzieva fc Herbstlll998l ). As shown in Fig. [TJ 
the PAH's have a strong enhancing effect on the formation of these molecules, especially on 
the more complex cyanopolyynes. Once again, there is little difference among the results for 
the three PAH models. In general, the dominant formation mechanism for each cyanopolyyne 
more complex than HC3N involves the next smaller member of the series. At 10^ yr, the 
HC5N, HC7N and HCgN abundances are enhanced by more than three orders of magnitude 
over their values without PAH's. The peak fractional abundances are indeed larger than the 
values observed in TM Cl (CP) of 2 x 10"^ 5 x 10"^° and 2.5 x 10"^° for HC5N, HC7N and 
HCgN respectively (see iMarkwick et al.ll2000l . and references therein). Similar enhancements 
are seen for other nitrogen-containing molecules such as CH2CHCN, CH3CN, CH3C3N and 
C3H2N+. 



4-2.5. sulfur-hearing species 

The elemental abundance of sulfur is smaller by more than 2 orders of magnitude in EAl 
compared with the two other sets of elemental abundances. As a consequence, most of the 
sulfur-bearing species have larger abundances with high-metal abundances at all times, with 
and without PAH's. In Fig. [HI we follow the temporal evolution of four species; OCS, H2S, CS 
and C3S. Three of these molecules tend to have enhanced peak abundances in the presence of 
PAH's, although the situation is more complex towards steady-state. The exceptional cases 
are CS for all elemental abundances and H2S for the low- metal case. The CS molecule is 
produced by the dissociative recombination of HCS"*" and HC2H"'" by electrons or dissociative 
neutralization of HCS"*" by PAH~. The larger ion HC2H"'" does not react with PAH~ to 
produce CS with our assumptions. As a consequence, the abundance of CS is smaller in the 
models including the PAH's. The late-time situation regarding H2S is more complex. The 
primary formation mechanism for H2S at most times is through the protonated precursor 
H3S"'", which is formed by a chain of reactions starting with S + H;^ and proceeding through 
HS+ and H2S+. Without PAH's, H2S is also produced by the charge exchange between H2S"'" 
and S. The ion H2S"'" is more abundant in the low-metal case without PAH's because it is 
produced by the reaction S+ + H2CO — > ^2^^ + CO where S''" is more abundant than the 
neutral form in the absence of PAH's. 
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5. Sensitivity to several PAH parameters 

There are some observational and theoretical indications that PAH's in dense clouds 
may be aggregated into clusters with radius up to 100 A. We have performed additional 
calculations using the three elemental abundances and PAH-clusters of radii 10 A and 100 
A with the reduced abundances discussed in Sect. |2i In our discussion of the results, we 
use the word PAH whatever the size. For the larger PAH's (10 and 100 A), we consider 
two expressions for the rate coefficients of the neutralization reactions between PAH^ and 
positive ions: one involving charge attraction, as expressed in equation (jlj), and the other 
a simple ballistic one, as expressed in equation The ballistic rate is considerably lower. 
For radii of 10 and 100 A and a temperature of 10 K, the ballistic rate is, respectively, about 
1000 X and 10 x smaller than the rate involving charge attraction. Another process with 
a rate coefficient that depends on the PAH size is the radiative attachment of electrons to 
neutral PAH's, which increases with increasing radius. In Fig. [9], we show the ionization 
fraction as a function of time with each set of elemental abundances and PAH's of 4, 10, 
and 100 A. As can be seen, the ionization fraction depends highly on the size of the PAH, 
and on whether or not charge attraction is considered. We do not show the dependence on 
the total abundance of C atoms incorporated into PAH's, but a larger abundance implies a 
more rapid decrease of the ionization fraction to a smaller value at late times. 

In the left panel of Fig. [HI neutralization proceeds with charge attraction for all three 
sizes of PAH's. For the two smaller sizes of PAH's, the main negatively-charged species in 
the model remains PAH~. For the PAH's of 100 A radius, the total abundance of PAH's is 
2 X 10~^^, much smaller than the abundance of electrons. The evolution of the ionization 
fraction is then closer to the case without PAH's, but the steady-state ionization fraction can 
still be somewhat lower. For the case with model EAl and the largest PAH's, the ionization 
fraction is the same with or without charge attraction and is not different from the model 
in which there are no PAH's at all. As a consequence, the 10 A and (for EA2 and EA3) the 
100 A cases represent intermediate stages between the models with 4 A PAH's and those 
without PAH's. The right-hand panel of Fig. M shows the results when positive ion-PAH~ 
neutralization occurs ballistically for the larger PAH's. In this case, the 10 A and 100 A 
cases give very similar results to the model without PAH's so that the effect of large PAH's 
at the chosen abundances is negligible in dense cloud chemistry. This effect occurs because 
the major depletion mechanism for PAH~ is photodetachment rather than recombination 
with positive ions, so that lowered neutralization rates have the sole effect of decoupling 
PAH^ from the ion-molecule chemistry. 

The modification of the ionization fraction has consequences for the abundances of 
species. In general, PAH's of 100 A (considering charge attraction or not) and of 10 A 
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without charge attraction lead to abundances similar to the models without PAH's. For 
most of the species, PAH's of 10 A with charge attraction give intermediate abundances 
between the case of very small 4 A PAH's and the case without PAH's. Many species such 
as C, Hj]", HCO"'" and OH show a difference smaller than a factor of two while using larger 
PAH's of 10 A instead of smaller PAH's of 4 A. Some species are largely underproduced 
using the elemental abundances EA2 and EA3 without PAHs. This is the case for C3H2, 
C3H4, C6H2, CH3C4H, CH3CHO and cyanoppolyynes. In the presence of 10 A PAH's, these 
species are almost as abundant as with smaller PAH's. 



6. Comparison with observations 



The objective of this section is to study the impact of PAH's on the comparison with 
observations in the two cold cores L134N (north peak) and TMC-1 (cyanopolyyne peak). 
These two clouds are known to show differences in the abun dances of C-rich mole cules that 
are sometimes explained by differences in age or C/0 ratio (IWakelam et al.ll2006l ). Here we 
explore the possibility that they are also explained by the amount and size of PAH's present. 



6.1. Method 



The observed abundances in L134N 



north peak) and TMC-1 (cyanopolyyne peak) are 



listed in Tables 3 and 4 of iGarrod et al.l ( 120071 ). In addition to these species, comprising 42 
molecules for L134N and 53 for TMC-1, we considered SiO, for which only upper limits o f 
3.6 X 10-^2 L134N and 2.4 x 10^^^ fo^ TMCl have been determined Jziurvs et al.lll989ah . 



In this comparison, we utilize the three different sets of elemental abundances listed in 
Table 1. Of all these elements, sulfur is the least constrained by abund ances of gas-phas e 
molecules gas-phase abundances so that the depleted fraction is uncertain (IRuffle et al.lll999l ). 
PAH size and abundance are also poorly determined parameters. Instead of fixing these 
parameters to previously used values, we try to give some trends by comparing observations 
with grids of models in which these three parameters are varied in addition to time. 

We computed 441 models by varying the abundance of elemental sulfur between 10^^ 
and 10~^, the abundance of PAH's between 10~^° and 10~^ and the size of PAH's between 4 
and 400 A. Note that the largest values of PAH sizes and abundances used here are not in 
agreement with the relation Nc = 460(a/10A)^ used in the first part of the article. Consid- 
ering the little information available on the relation of size to number of carbons for PAHs, 
in particular for PAH clusters, we consider the PAH size and abundance as independent. 
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The agreement between present models and observations quantified using statistical com- 
parisons, minimization for instance, is not good enough to give reasonable constraints on 
model parameters. For this reason, we did not use such methods but instead for each model 
and each time, we count the number of species with fractional abundance X reproduced by 
the model within one order of magnitude (Xobs/10 < Xmodei < Xobs x 10) where Xobs is the 
observed abundance and Xmodei the computed abundance. For those species for which only 
upper limits to their abunda nces are known a nd for the case of atomic carbon, for which 
only a lower limit is known (jStark et al.lll996l ). we check that the predicted abundance is 
in agreement with the observed limit when this is the case, we consider that the species is 
reproduced. 



6.2. Results 

Figures [TD] and [TT] show the maximum percentage of species reproduced in each cloud 
as a function of the parameters chosen. Column (a) shows the level of agreement as a 
function of PAH size and abundance, with each point optimized for the best elemental 
abundance of sulfur and the best time. Column (b) shows the level of agreement as a 
function of the elemental abundance of sulfur and the time optimized at each point for the 
best PAH parameters. Finally, column (c) shows the agreement as a function of the elemental 
abundance of sulfur and the time for models without PAH's. These results are shown for 
each set of elemental abundances (1: EAl, 2: EA2 and 3: EA3). Best agreement, with or 
without PAH's, is most often achieved at so-called early times of 10^~^ yr. Moreover, high 
levels of agreement cannot be achieved with high abundances of PAH clusters, as shown in 
the upper-right corners in Column 1 for both figures. 

The maximum percentage of molecules reproduced depends on all the parameters, and 
several solutions exist. PAH's do not improve significantly the best agreement for L134N 
(already above 70% without PAH's) but they do expand the range of parameters (time and 
elemental abundance of sulfur) for which we have an agreement larger than 70%, especially 
if elemental abundances 1 or 2 are utilized. Note particularly the expanded range of agree- 
ment with the use of elemental abundances 2 (the high-metal ones). The optimum PAH 
parameters corresponding to the yellow zones in Column b do not depend on the optimum 
sulfur abundance. As an example, high agreement is obtained with 3 x 10~^° abundance of 
4A PAH's both for low (3 x 10~*) and high (1 x 10~^) elemental abundance of sulfur. The 
chemistry in TMC-1 on the contrary is better reproduced if PAH's are present, as can be seen 
immediately by comparison of columns b and c in Figure [TTl Without PAH's, agreement 
rarely extends to more than 60%, while with PAH's, agreement in excess of 70% exists for all 
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three sets of elemental abundances. For elemental abundances 1, the PAH parameters are 
relatively unconstrained; while for the high-metal abundances (EA2), high abundances of 
small PAH's are needed. The situation is less clear for EA3. For both clouds, we find isolated 
areas of agreement larger than 70% at very early times, < 1000 yr, for large abundances of 
sulfur (> 8 X 10"^) and big (> 25A) PAH's (> 10"^) if EA2 and/or EA3 is utilized. Such 
agreement is of mathematical interest only. 



7. Conclusion 

Most chemical studies of the dense interstellar medium rely purely on gas-phase pro- 
cesses. A few contain, in addition, a varying amount of surface chemistry on interstellar 
grains but the inclusion of PAH's, which occupy a middle ground in size between small 
gas-phase species and dust particles, is quite rare. Part of the reason for this neglect is 
that to include PAH's properly, one should view them as individual species. But the ac- 
tual individual species that make up the PAH content of interstellar clouds are not well 
constrained from spectral observations, and even their chemistry has not been studied in 
great detail. So, inclusion of PAH's means that one must treat them as small grains, and 
follow just their charge state. In this study, as in previous ones, we lump all PAH's into 
just three groups, depending on whether they are neutral, positively-charged, or negatively- 
charged. Positively-charged PAH's were initially assumed to be unimportant in cold dense 
sources, an assumption that was confirmed upon the request of the referee. They were 
not included in the chemical network used for the results presented here. By including an 
assortment of chemical processes involving negatively-charged and neutral PAH's into our 
standard gas-phase network of reactions, we are able to determine the effect of PAH's on 
the time-dependent gas-phase chemistry of cold and dense interstellar cloud cores. The 
gas-phase model employed is known as osu_01_2007, and can be downloaded from the 



URL |http://www.physics.o hio-state.edu/^eric/research.html To this model, we add elec- 



tron attachment to neutral PAH's, neutralization between negatively- charged PAH's and 
positively-charged smaller ions, and photodetachment of negatively-charged PAH's by exter- 
nal photons. Our standard model contains a fractional abundance of 3.07 x 10"'' PAH's of 
size 4 A, corresponding to 30 carbon atoms. In addition to this model, we vary the size and 
abundance of the PAH's to study the impact of these parameters on the chemistry. We also 
use three different sets of elemental abundances with different abundances of heavy atoms. 

In general, the dominant effect in our calculations is that electrons are replaced by PAH~ 
as the dominant negative charge carrier, an effect which weakens and eventually vanishes as 
the size of the PAH's, due probably to clustering, and their fractional abundance decrease. 
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When the abundance of PAH~ exceeds that of electrons, the overall ionization fraction is 
reduced because atomic positive ions are far more rapidly neutralized by reactions with PAH~ 
than by radiative recombination with electrons. These changes lead to many secondary 
effects. Our calculations show that many of the molecular species have abundances that 
are strongly affected by the presence of PAH's. In particular, molecular ions tend to have 
increased abundances at most times and large unsaturated neutrals, such as the polyynes, 
as well as saturated species such as methanol, show very large increases in abundance. The 
sensitivity of our results to variations in the size of PAH's depends upon assumptions made in 
the rate of neutralization of PAH~ and positive ions. The differences among results obtained 
with different elemental abundances are much smaller in models containing PAH's than in 
models without them. 

By comparing observations in the two dense cloud cores L134N and TMC-1 with grids 
of models, we showed only that the chemistry in such sources is not consistent with large 
amounts of large PAH clusters. We are not able to give more severe unambiguous constraints 
on the PAH parameters, since high levels of agreement can exist for large parameter varia- 
tions. The situation is similar for the elemental abundance of sulfur, where values over several 
orders of magnitude are consistent with maximum agreement with observations. As regards 
time, both the standard early-time and very early times can lead to optimum agreement, 
although the latter must be rejected on astronomical grounds. The situations for L134N and 
for TMC-1 show differences. For L134N, the best agreement is not improved when PAH's 
are included, although the range of optimum sulfur elemental abundances and times are 
extended, most significantly so for the high- metal abundances (EA 2). For TMC-1, on the 
other hand, the best level of agreement is significantly better when PAH's are included, using 
any of the three elemental abundances. This difference may indicate a difference in PAH 
properties between the two regions. In particular, for EA 2, optimum agreement in TMC-1 
is obtained with smaller and more numerous PAH's, while for L134N, such parameters make 
matters worse. Finally, unlike the case of gas-phase models without PAH's, models including 
PAH's are able to accommodate high-metal elemental abundances, similar to those observed 
in the diffuse medium, if the effects of PAH's are included and optimized. 

V. W. acknowledges partial financial support from the CNRS/INSU program PCMI. E. 
H. thanks the National Science Foundation for its partial support of this work. 
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Species 


EAl 


1 


EA2 


EA3 


Reference for EA2 


He 


1.40( 


-1) 


9.00(-2) 


9.00(-2) 


see text 


N 


2.14( 


-5) 


7.60(-5) 


7.60(-5) 


Cardelli et al. ri993) 





1.76( 


-4) 


2.56(-4)^ 


2.56(-4)^ 


Mever et al. fl998) 


c+ 


7.30( 


-5) 


1.20(-4)2 


1.20(-4)2 


Cardelli et al. ri993) 


s+ 


8.00( 


-8) 


1.50(-5) 


1.50(-5) 


Federman et al. fl993)^ 


Si+ 


8.00( 


-9) 


1.70(-6) 


0.00 


Cardelli et al. ri994) 


Fe+ 


3.00( 


-9) 


2.00(-7) 


1.50(-8) 


Savasre et al. ri992) 


Na+ 


2.00( 


-9) 


2.00(-7) 


2.00(-7) 


Savaee & Sembach f 19961 


Mg+ 


7.00( 


-9) 


2.40(-6) 


0.00 


Savaee et al. fl992) 


C1+ 


1.00( 


-9) 


1.8(-7) 


1.80(-8) 


Savaee & Sembach fl996')'^ 


P+ 


2.00( 


-10) 


1.17(-7) 


1.17(-7) 


Savaee & Sembach fl996l^ 


F+ 


6.68( 


-9)^ 


1.8(-8) 


6.68(-9)5 


Federman et al. (2005) 



a(-b) refers to a x 10~^ 

^ow-Mctal abundances from ICraedel et all |l982h : iMortonI |l974l ) 

^We considered an additional depletion compared to the observed value to take into account the oxygen and 
carbon depleted on grain mantles (see text). 

^The sulfur abundance found by the authors is larger than what is usually assumed as the cosmic abundance. 
We use the cosmic abundance of 1.5 x 10^^. 



Abundances have been recalculated from the observations of iMortonI (|l974[ ). 
5 Depleted value from lNeufeld et all |2005h . 
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Fig. 1. — Ionization fraction as a function of time for the three elemental abundances (solid 
line: EAl; dashed line: EA2; dashed-dotted line: EA3) without PAH's (grey lines; left panel) 
and with PAH's (black lines; right panel). 
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Fig. 2. — Electron, PAH and PAH abundances plotted vs. time for the model EA2/PAH. 
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Fig. 3. — C and abundances as a function of time for the three elemental abundances 
(solid line: EAl; dashed line: EA2; dashed-dotted line: EA3) without PAH's (grey lines) 
and with PAH's (black lines). 
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Fig. 4. — B.^ , HCO^, H3O"'", OH, H2O and O2 abundances as a function of time for the 
three elemental abundances (solid line: EAl; dashed line: EA2; dashed-dotted line: EA3) 
without PAH's (grey lines) and with PAH's (black lines). 
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Fig. 5. — C-C3H2 (cyclopropenylidene), C3H4 (methyl acetytene), C6H2 (triacetylene) and 
CH3C4H (methyl diacetylene) abundances as a function of time for EAl (solid lines), EA2 
(dashed lines), and EA3 (dashed-dotted lines) without PAH's (grey lines) and with PAH's 
(black lines). For the three last molecules, the abundance produced by model EA2 and EA3 
without PAH's is smaller than 10~^^. 
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Fig. 6. — The abundances of formic acid (HCOOH), ketene (CH2CO), methanol (CH3OH), 
and acetaldehyde (CH3CHO) are plotted as functions of time for three elemental abundances 
(solid line: EAl; dashed line: EA2; dashed-dotted line: EA3) without PAH's (grey lines) and 
with PAH's (black lines). The acetaldehyde abundance in model EA2 and EA3 is smaller 
than 10"^^ 
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Fig. 7. — Cyanopolyyne abundances as a function of time for the three elemental abundances 
(solid line: EAl; dashed line: EA2; dashed-dotted line: EA3) without PAH's (grey lines) 
and with PAH's (black lines). The HC7N and HCgN abundances with models EA2 and EA3 
are too small to be seen. 
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Fig. 8. — The abundances of OCS, H2S, CS and C3S as functions of time for the three 
elemental abundances (solid line: EAl; dashed line: EA2; dashed-dotted line: EA3) without 
PAH's (grey lines) and with PAH's (black lines). 
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Fig. 9. — Ionization fraction as a function of time for small and large PAH's using the three 
elemental abundances. For the two larger PAH's, we use expressions for the rate coefficients 
of PAH~-X"'" neutralization reactions that include (left column) and exclude charge attraction 
(right column). 



- 32 - 



I: 



■ggg 



1 


I 


1 


I 


1 


1 


] 


■ 


■ 


B 


m 


m 








ir 






m 






'I 


■ 




m 


m 


m 




w 














■ i. . 















_ifl _s _3 _7 





1 


1 


t 


i 


a 
1 


1 


1 

li 




I 






■ 


M 




■ 


'•71 




11 


■ 


• 


m 


■ 


■ 




■1 






* 


4 






^ 10. 


« 















3-a 



5; 



10 



] I I I t 1 

■ t m m m % 

m k If H b « 



* m 



mm*' 



□ 



-a 

LT 

? -4 



^ I 



■|ci^(i::nB> 

2-b 



7; « 




lcig(!lm£i} 



□ 

■s 

□ 



■b 

1^ 



-D 



I — 



Ing(lh^) 



Fig. 10. — Level of agreement (violet: less than 50%; blue, between 50 and 60%; green, 
between 60 and 70%; and yellow, between 70 and 80%) between models and observations in 
L134N. Each row represents the results for a set of elemental abundances (1: EAl, 2: EA2 
and 3: EA3). Column a represents the level of agreement as a function of PAH parameters 
with optimized sulfur and time at each point, best agreements for sulfur elemental abundance 
and time as a function of PAH parameters. Column b represents the level of agreement as 
a function of the sulfur elemental abundance and time with optimized PAH parameters at 
each point. Column c represents the agreement without PAH's. 
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Fig. 11. — Level of agreement (violet: less than 50%; blue, between 50 and 60%; green, 
between 60 and 70%; and yellow, between 70 and 80%) between models and observations in 
TMC-1. Each row represents the results for a set of elemental abundances (1: EAl, 2: EA2 
and 3: EA3). Column a represents the level of agreement as a function of PAH parameters 
with optimized sulfur and time at each point, best agreements for sulfur elemental abundance 
and time as a function of PAH parameters. Column b represents the level of agreement as 
a function of the sulfur elemental abundance and time with optimized PAH parameters at 
each point. Column c represents the agreement without PAH's. 



